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Abstract 

In this paper, the filter-and-forward (FF) relay design for orthogonal frequency-division multiplexing (OFDM) 
transmission systems is considered to improve the system performance over simple amplify-and-forward (AF) relaying. 
Unlike conventional OFDM relays performing OFDM demodulation and remodulating, to reduce processing complex- 
ity, the proposed FF relay directly filters the incoming signal with a finite impulse response (FIR) and forwards the 
filtered signal to the destination. Two design criteria are considered to optimize the relay filter for given source power 
allocation. One is to minimize the relay transmit power subject to per-subcarrier signal-to-noise ratio (SNR) constraints 
and the other is to maximize the worst subcarrier channel SNR subject to the total relay transmit power constraint. It is 
shown that the two problems reduce to semi-definite programming (SDP) problems. Furthermore, the problem of joint 
source power allocation and relay filter design is considered for the second criterion, and a convergent iterative algo- 
rithm based on alternating optimization is proposed for the joint design. Numerical results show that the proposed FF 
relay significantly outperforms simple AF relays under both criteria. Thus, the proposed FF relay provides a practical 
alternative to the AF relaying scheme for OFDM transmission. 
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I. Introduction 

Recently, relay networks have drawn extensive interest from the research community because they play 
an important role in enlarging the network coverage and improving the system performance in current and 
future wireless networks. Indeed, LTE-Advanced employs relays for coverage extension and performance 
improvement [1]. There are several well-known relaying schemes such as AF, decode-and-forward (DF), 
and compress-and-forward (CF) [2-4]. Among the relaying schemes, the AF scheme is the simplest and 
is suitable for cheap relay deployment under transparent* operation [1]. Recently, there have been some 
efforts to extend this simple AF scheme to a linear filtering relaying scheme, i.e., an FF scheme, to obtain 
better performance than the AF scheme while keeping the benefit of low computational complexity of the AF 
scheme [5-10]. It has been shown that the FF scheme can outperform the AF scheme considerably. However, 
most of the previous works on the FF relay have been done for single-carrier transmission, whereas most of 
the current wireless standards adopt OFDM transmission. Thus, in this paper, we propose direct FF relaying 
for OFDM transmission instead of using conventional OFDM relays which OFDM-demodulate the incoming 
signal, amplify or decode the demodulated signal, OFDM-remodulate the processed signal and transmit the 
remodulated OFDM signal to the destination [11-14]. In the proposed scheme, the incoming signal to the 
relay is FIR filtered at the chip rate of the OFDM modulations and the filtered signal is directly forwarded 
to the destination. In this way, the necessity of OFDM processing at the relay is eliminated, but the overall 
performance can still be improved over the AF scheme by a properly designed relay filter. 
A. Our Approach and Contributions 

There are several design criteria for the FF relay design for OFDM systems: rate maximization, minimizing 

power consumption, or overall quality-of-service (QoS) improvement. In this paper, we mainly consider 

power consumption minimization and QoS improvement under the scenario of single-input single-output 

(SISO) communication as the first step in this research direction.^ First, we consider the case that the source 

power allocation is given and the FF relay filter is optimized. In this case, we consider two design criteria. 

One is to minimize the relay transmit power subject to per-subcarrier SNR constraints, and the other is to 

maximize the worst subcarrier channel SNR subject to a total relay transmit power constraint. The second 

criterion is particularly interesting from the perspective of overall QoS over subcarrier channels. Under the 

assumption of known channel state between the source and the relay and known channel statistic between 

* Transparent operation means that the destination node does not know the existence of the relay node, and this operation is suitable 
for cheap AF relays [1]. 

^The multiple-input multiple-output (MIMO) communication case is beyond the scope of this paper and will be studied in a future 
work. For FF relay design for rate maximization, see Section V and the appendix. 
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the relay and the destination, which is suitable for cheap transparent relay operation, we show, by exploiting 
the eigen-property of circulant matrices and the structure of Toeplitz filtering matrices, that each problem can 
be reformulated as an SDP problem based on a semi-definite relaxation approach. Furthermore, for the first 
criterion, we show that the solution of the relaxed SDP problem is the same as that of the original problem 
under certain mild conditions. Next, we consider the case that both source power allocation and relay filter 
are optimized to maximize the worst subcarrier channel SNR under power constraints both at the source and 
the relay. For this joint optimization, we propose an iterative algorithm based on alternating optimization 
between linear programming for source power allocation and a SDP programming for the relay filter design, 
which converges to a local optimum at least. Numerical results show that the proposed FF relay significantly 
outperforms simple AF relays under both criteria. Thus, the proposed FF relay provides a practical alternative 
with low complexity to the AF relaying scheme for OFDM transmission. 
B. Notation and Organization 

In this paper, we will make use of standard notational conventions. Vectors and matrices are written in 
boldface with matrices in capitals. All vectors are column vectors. For a scalar a, a* denotes its complex 
conjugate. For a matrix A, A*, A T , A^, and tr(A) indicate the complex conjugate, transpose, conjugate 
transpose, and trace of A, respectively. A y and A y mean that A is positive semi-definite and that A is 
strictly positive definite, respectively. For two matrices A and B, A y B means that A — B y 0. I n stands 
for the identity matrix of size n (the subscript is omitted when unnecessary), and rnxn denotes a m x n matrix 
of all zero elements. The notation Toeplitz(f T , N) indicates a N x (N + Lf — 1) Toeplitz matrix with 
N rows and [f T , 0, • • • , 0] as its first row vector, where f T is a row vector of size Lt, and diag(<ii, • • • , d n ) 
means a diagonal matrix with diagonal elements d\, ■ ■ ■ , d n . The notation x ~ CJ\f{n, S) means that x 
is complex circularly-symmetric Gaussian distributed with mean vector /i, and covariance matrix X. 
denotes the expectation, j = \/— 1. 

The remainder of this paper is organized as follows. The system model is described in Section II. In 
Section III, the FF relay design problems are formulated and solved by using convex optimization theory. 
The performance of the proposed design methods is investigated in Section IV, followed by the conclusion 
and discussion in Section V. 

II. System Model 

We consider a full-duplex relay network composed of a source node (basestation), a relay node and a 
destination node (terminal station), as shown in Figures 1 and 2, where the source employs OFDM modulation 
with N subcarriers, and each link performs SISO communication. We consider the case that the direct link 
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between the source and the destination is weak. Thus, for simplicity, we assume that there is no direct 
link between the source and the destination and that both the source-to-relay (SR) link and the relay-to- 
destination (RD) link are inter-symbol interference (ISI) channels modelled as FIR filters. We assume that 
the relay is an FF relay, i.e., the relay performs FIR filtering on the incoming signal at the chip rate of the 
OFDM modulation and transmits the filtered output immediately to the destination. Thus, the FF relay can be 
regarded as an additional ISI channel between the source and the destination. We assume that the introduction 
of FIR filtering at the relay does not make the length of the overall FIR channel between the source and the 
destination larger than that of the OFDM cyclic prefix. (We will see later in Section IV that only a few taps for 
the relay FIR filter are enough to obtain most of the gain of the FF relay over the AF relay.) Since our focus 
of using the FF relay in this paper is the transparent relay operation, we assume that the SR channel state is 
known to the relay and that the RD channel state is unknown to the relay but the RD channel distribution is 
known to the relay. * Such an assumption is reasonable for the transparent relay operation since the destination 
node does not know the existence of the relays; thus, the destination node does not feedback information to 
the relay node directly. 

Specifically, at the source, the length N data vector of OFDM symbols is given by s := [s[0], s[l], • • • , 
s[N — l]] T , where each data symbol is assumed to be a zero-mean independent complex Gaussian random 
variable, i.e., s[k] ~ CM(0, P s ,k) f° r k = 0, 1, • • • , N — 1. The time-domain vector x s after normalized 
inverse discrete Fourier transform (IDFT) at the source is given by 



x,[0] 
x s [l] 


1 


1 
1 


1 


1 

JV-1 
N 




s[0] 
s[l] 


x a [N-l] _ 




1 








s[N - 1] 


= :x s 






=:W N 






= S 



(1) 



where oj^ = e? ' n . Let wTi denote the k-th row of the normalized IDFT matrix Wjv for k = 1, ■ ■ ■ , N. 
Then, x s can be written as x s = [wqS, w^s, • • • , ■wJ f _ 1 s] T and the covariance matrix S Xs of x s is given 

*The channel state information for the SR channel can be obtained at the relay by channel estimation based on possible training 
or a preamble signal embedded in the transmitted signal from the source. The channel statistic information for the RD link can be 
obtained based on the size of the relay cell or some initial measurement when the relay is deployed. It will be shown later that the 
proposed FF relay design is robust against channel statistic mismatch for the RD channel. 
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by 



E{x,x, 



J_ T-^N-l p 



™N-2Ps 



N l^k=0 r s,k 



(2) 



where p s = [P a , P s l , • • • , P s , N -i] , since 



E{w/ss"w*} 



1 



if i > j, 



(3) 



if i < j. 



The vector x s is attached by a cyclic prefix with length Lcp, i.e., 



x,n 



n = 0,1, ••• ,JV- 1, 



n 



1,-2,- •• 



(4) 



and the cyclic prefix attached sequence x s [n] is transmitted from the source to the relay through the SR 
channel. Then, the received baseband signal at the relay is given by 

L f -1 



y 



rN = ^2 h% s [n -l]+ n r [n], 



(5) 



1=0 



where f = [/ , A, • • • , /l 7 -i] T is the channel tap coefficient vector of the SR channel known to the relay, 
Lf is the length of the SR FIR channel, and n r [n] is the additive white Gaussian noise at the relay with 
n r [n] ~ C/V(0,of )■ At the relay, the received signal y r [n] is FIR filtered at the chip rate of the OFDM 
transmission and then transmitted immediately to the destination. Thus, the output signal at the relay at 
(chip) time n is given by 

L r -i 

Vt[n] = ^2 n Vr[n - I], (6) 

1=0 

where r = [tq, r\, ■ ■ ■ , r£ r _i] T is the FIR filter coefficient vector at the relay and L r is the order of the FIR 
filter. Note that, when L r = 1, the FF relay simply reduces to the AF relay. However, when L r > 1, the FF 
relay is an extension of the AF relay with some amount of digital processing. Finally, the signal transmitted 
by the relay goes through the RD FIR channel to the destination; thus, the received signal at the destination 
is given by 

Lg-l 

Vd[n] = ^2 9i Vt[n + n d [n], (7) 

1=0 
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where g = [go, gi, ■ " > 5L g -i] T is the FIR channel tap coefficient vector for the RD channel, L g is the order 
of the RD FIR channel, and n d [n] is zero-mean white Gaussian noise with variance a\ at the destination. 
Here, we assume that the channel tap coefficient gi, I = 0, 1, • • • , L g — 1, is independent and identically 
distributed (i.i.d.) according to g\ 1 ~ ' CN(0, a g ), i.e., each tap is independently Rayleigh faded, and that the 
realization {gi, I = 0, 1, • • • , L g — 1} is not known to the relay but its distribution is known to the relay. By 
stacking the output symbols at the relay and the received symbols at the destination, we have the following 
vectors for the transmitted signal at the relay and the cyclic prefix portion removed received signal vector at 
the destination, respectively: 

y t = RFx s + Rn r and y d = GRFx s + GRn r + n d , (8) 

where 



yd 


= [y d [N-l],y d [N-2], 


••• ,y d .M] T , 


(9) 


yt 


= [y t [N-l},y t [N-2],- 


•• ,y t [Q},y t [-l},--- ,y t [-L g + l]] T , 


(10) 


x. s 


= [x a [N-l],x,[N-2], 


■■■ ,x s [Q],x s [-l], ■ • • ,x s [-L g - L r - L f + 3]] T , 


(11) 


n, 


= [n r [N-l],n r [N-2], 


••• ,n r [0],n r [-l] • • • ,n r [-L g - L r + 2]] T , 


(12) 


n d 


= [n d [N-l],n d [N-2] 


, n dM\ T , 


(13) 


G 


= Toeplitz(g T , N), 




(14) 


R 


= Toeplitz(r T , N + L g 


-i), 


(15) 


F 


= Toeplitz(f T ,iV + L g 


+ L r -2). 


(16) 



Under the assumption that Lqp > L g + L r + Lf — 3, the DFT of the cyclic prefix portion removed received 
vector of size N at the destination is given by 

y d = W#GRFx s + W#GRn r + W#n d , (17) 
= W§U c W N s + W$GRn r + W%n d , (18) 
= Ds + Wf GRn r + n d , (19) 

where is the normalized DFT matrix of size N, H c is a N x N circulant matrix generated from the 
overall Toeplitz filtering matrix GRF from the source to the destination, and D = diag(cio, • • • , djv-i) = 
W^H c Wtv is the eigen-decomposition of H c . 
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A. Manipulation for quadratic forms 

The received signal form (19) is standard in OFDM transmission, but the form cannot be used directly for 
relay filter optimization in the next section. Thus, in this subsection, we derive an explicit expression for the 
received signal y d [k], k = 0, 1, • • • , N — 1, which facilitates optimization formulation in the next section, 
based on the following property of circulant matrices [15]. 

Lemma 1: [15] Let C be an N x N circulant matrix with the first row [c(0), c(l), • • • , c(N — 1)] . Then, 
the eigenvalues of C are given by 

N-l 

A fc = ^c( n )^", fc = 0,V ,N-1, 

n=0 

with the corresponding right eigenvectors 

£ fc = -jL[l,w^o#V- ,oj- {N - 1)k ] T , fc = 0,l,.-. ,N-1. 
V iV 



(21) 



By Lemma 1, to derive the diagonal elements of D in (19), we need to know only the first row of H c in (18). 
Let the first row of G be g T . Then, g T is a 1 x (N + L g — 1) row vector given by 

g T = [g T , (),••• ,0], (20) 

and the first row of GRF is given by g T RF. Since H c is generated by truncating out the elements of GRF 
outside the first N x N positions and by moving the lower (L g + L r + Lj — 3) x (L g + L r + Lf — 3) elements 
of the truncated part to the lower left of the untruncated N x N matrix, the first row h^T of H c is simply the 
first N elements of the first row g r RF of GRF, i.e., h^T = g r RFT, where T is a truncation matrix for 
truncating out the elements of g T RF except the first N elements, given by 

0(Lg+L r +L f -3)xN 

Now, the diagonal elements of D can be obtained by Lemma 1 and are given by 

[d , ■ ■ ■ ,ci7V-i] T = ViVW$(g T RFT) T , (22) 

where V^W^ is the DFT matrix of size 7Y. Finally, the received signal in the k-th subcarrier at the destina- 
tion is expressed as 

y d [k] = ViVwf T T F T R T g S [fc] + wfGRn, + wf (23) 
where is the (k + l)-th row of Wjy. Thus, the signal and noise parts of yd[k] are given by 

Vd,s[k] = VNw% T T F T R T gs[A;] and y d>N [k] = wfGRn r + wf n d , (24) 
respectively, for k = 0, 1, ■ ■ ■ , N — 1. 
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III. FlLTER-AND-FORWARD RELAY DESIGN CRITERIA AND OPTIMIZATION 
In this section, we consider several meaningful FF relay design problems for the relay network employing 
OFDM transmission described in Section II. First, we consider the case that the source power allocation 
{Ps,k} is g lVen > and the FF relay filter is optimized. In this case, we consider two design criteria. The 
first criterion is to minimize the transmit power of the FF relay subject to a SNR constraint for each OFDM 
subcarrier channel, and the second criterion is to maximize the worst subcarrier SNR subject to a relay 
power constraint to improve overall QoS for subcarrier channels. Here, we shall show that each problem can 
be reformulated as a semi-definite programming (SDP) problem based on a semi-definite relaxation (SDP) 
approach. That is, the original non-convex FF relay design problems are approximated by convex SDP 
problems. For the first criterion, we shall show that the solution of the relaxed SDP problem is the same as 
that of the original problem under certain mild conditions. Next, we consider the case that both source power 
allocation and relay filter are optimized to maximize the worst subcarrier SNR under power constraints both 
at the source and the relay. For this joint optimization, we propose an iterative algorithm based on alternating 
optimization between linear programming for source power allocation and SDP programming for the relay 
filter design, which converges to a local optimum at least. 

A. FF relay optimization for given source power allocation 

A. 1 FF relay transmit power minimization under per-subcarrier SNR constraints 

We first consider the problem of designing the FF relay tap coefficient vector r = [ro,rx, ■ ■ ■ , rL r _i] T 
to minimize the relay transmit power subject to an SNR constraint per OFDM subcarrier. This problem is 
formulated as follows: 

Problem 1: For given source power allocation {P S: k, k = 0, 1, ■ ■ ■ , N — 1}, SR channel f, RD channel 
static information {a^,L g ), FF relay filter order L r , and a set {7^, k € 1} of minimum desired SNR values 
for subcanier channels X, 

min P r subject to (s.t.) SNR fc > 7*., V k G 1 C {0, 1, • • • , TV - 1} (25) 

r 

where P r is the relay transmit power and SNR& is the SNR of the fc-th subcarrier channel. 

To solve Problem 1, we need to express each term in the problem as a function of the design parameter r. First, 
let us derive the SNR on the /c-th subcarrier channel in the received signal (23) at the destination. Note that 
the signal and noise parts in (24) are represented in terms of the relay filtering matrix R. The representation 
of SNR in terms of R is redundant since r is embedded in R. Thus, we need reparameterization of SNR in 
terms of r, and this is done based on (24) by exploiting the Toeplitz structure of R as follows. Using (24), we 
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n\y d ,sM 2 } ■ 



(c 

(/) 



iVE{wfT T F T R r g|s[/c]| 2 g^R*F*T*w A .}, 
iVwfT T F T R T IE{|s[A;]| 2 gg j:/ }R*F*T*w fc , 
iYP S)fe tr(wfT T F T R T E{gg ff }R*F*T* WS; ) ; 
AfP s , feC T|tr(wfT T F r R T I L9 I Lfl R*F*T*w fc ) 
iVP s , fcC T 5 2 tr(F*T*w fc wfT T F T R^R2 s ), 



iVP Sifc a 2 tr(R2 9 K fc Ri 9 ), 



r 

H 



vec(RiJ K fc vec(Ri o 



iVP Sifc a 2 [vec(R T )] H K fc vec(R 
iVP a>fc ^ r ff El K fc Efr, 



where 



Ii 9 °L 9 x(Ar-l) 
0(JV-l)xL 9 0(jy-l)x(iV-l) 





r l 9 




; R = 




; K fc 




Rjv-i 





Rl s = Il 9 R is a matrix composed of the first L g rows of R; and 

El = [ Ii r , irX (7v+L 9 -2) > 0-LrXl, Ii r ) 0L r x(Ar+L 9 -3)) " ' " > 0i r x (N+L g -2) > I.Lr 



Af+L 9 +L r -2 columns N +L„+L r —2 columns 



Af+L„+L r -2 columns 



(26) 



IiJK t ; K fc = Ir <g)K fc ; 



(27) 



Here, (a) is due to the assumption of independence of the signal and the RD channel coefficients; (b) is due to 
the assumption^ of g x ' CM(0, cr 2 ) (see (20)) ; (c) and (d) are due to tr(ABC) = tr(CAB); (e) is due to 



H 



tr(R£ KfcRf ) = vec(R£ 9 ) K fc vec(R| /g ); and (f) is obtained because R = Toeplitz(i~ z , N + L g - 1) 

and thus vec(R r ) = E^r. The key point of the derivation of (26) is that the received signal power at the fc-th 

subcarrier channel is represented as a quadratic form of the design variable r. Next, consider the received 

noise power for the fc-th subcarrier channel. Using similar techniques to those used to obtain (26), we can 

§ The i.i.d. assumption for gi is not necessary and the assumption of independence is enough. Under the independence assumption, 
however, the formulation becomes more complicated. 
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express the received noise power based on the noise part in (24) as 



n\yd,NM 2 } = (T r 2 tr(R fl G%wfGR) + a 



di 



a z r [vec(R)] H M fc vec(R) + o%, 



a 2 r r H E 2 M k E$r + aj, 



(28) 



where M k = I N+L +Lr -2 <8> M k and E 2 is given by 



El 



1 e( e 
T 



3 T 

-N+L a 



-i o* 



T 
T 



J 



-iV+L„-l 



(29) 



Here, ej be the i-th row of Ijv+L g -i and the size of each T in (29) is 1 x (N + L g — 1). (It is easy to 
verify that vec(R) = E^r due to the Toeplitz structure of R.) Based on (26) and (28), the SNR of the fe-th 
subcarrier channel is now expressed as 



SNR, 



ATP^r^EiKfcEfr 



(30) 



E 2 M fc E^r + a\ 

Finally, consider the relay transmit power. Using (8), we obtain the relay transmit power in a similar way as 



E{tr(y t yf)} 



tr(RFE{i s if } F ff R H ) + tr^RR^), 
tr(R(F£ is F" + a r 2 I)R H ), 



= : n 



[vec(R H )] flvec(R 



r T EiflEf r* = r^EilTEf r, 



(31) 



where n = In + l 9 -i ® n, and ~E^ a is obtained similarly to (2) based on (3). The last equality holds since 
the power is a real-valued quantity. 

Now, based on (26), (28) and (31), Problem 1 can be restated as follows: 



Problem 2: For the given parameters, 



min r^EilTEf r s.t. 1 - > lk , kel. 



(32) 



March 4, 2013 



DRAFT 



SUBMITTED TO IEEE TRANS. ON WIRELESS COMMUNICATIONS, March 4, 2013 11 

In general, Problem 2 is not a convex problem. However, the problem can still be solved efficiently by using 
convex optimization techniques. Let 1Z := rr H . Then, by using tr(ABC) = tr(BCA) and relaxing the rank 
one constraint for TZ, Problem 2 can be reformulated as follows: 

Problem 3: For the given parameters, 

min tr($ P lZ) 
n 

s.t. tr ([*s(fc) - lk 3> N {k)]1Z) > aj lk , k€T, 

TZhO (33) 
where * P = EjH'Ef , *<?(*;) = TVP^EiKfcEf , $ N (k) = a r 2 E 2 M fe Ef . 

Note that by relaxing the rank one constraint for TZ, Problem 2 is converted to Problem 3, which is a semi- 
definite program (SDP) [16] and it can be solved efficiently by using the standard interior point method for 
convex optimization [16], [17]. With an additional constraint rank (TZ) = 1, Problem 3 is equivalent to the 
original Problem 2. That is, if the optimal solution to Problem 3 has rank one, then it is also the optimal 
solution to Problem 2. However, there is no guarantee that an algorithm for solving Problem 3 yields the 
desired rank one solution. In such a case, randomization techniques [18] can be used to obtain a rank-one 
solution r from TZ. However, for this specific problem related to the transparent FF relay design, we provide 
a stronger result, stated in the following theorem. 

Theorem 1: If all the desired SNR constraints except one are satisfied with strict inequality, then the non- 
trivial optimal solution of Problem 3, which is a relaxed version of Problem 2, always has rank one. 
Proof : See the appendix. 

Note that the condition in Theorem 1 is mild and is satisfied in many cases. Thus, solving Problem 3 directly 
yields the solution to the original power minimization problem under subcarrier SNR constraints in many 
cases. 

A.2 Worst subcarrier SNR maximization under a relay transmit power constraint 

Now, we consider another useful FF relay design problem for given source power allocation. The second 
problem is to maximize the SNR of the worst subcarrier channel under a total relay transmit power constraint. 
This problem is formulated as follows: 

Problem 4: For given source power allocation {P S: k, k = 0, 1, ■ ■ ■ , N — 1}, SR channel f, RD channel 
static information (cr 2 , L g ), FF relay filter order L r , and maximum available relay transmit power P r , m ax» 

max min SNRjt s.t. P r < P rmax . (34) 

r ke{0,- ,N-1} 
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By using (30) and (31), Problem 4 can be expressed as 

NP S fecr^EiKi-Ef r 
max min - g ' 9 k 1 s.t. r^EjlTEf r < P r , m3X . (35) 

r fce{o,-,Ar-i}ofr ff E 2 M fc E^r + aj 

By introducing a slack variable r, the above problem can be rewritten as 

max r 

r 

NP S fefT^EiKjuEfr 

s.t. 1 - > r , fc = 0, 1, . . . , JV - 1 

a r 2 r"E 2 M fc E^r + a\ ~ ' 

r^EilTEf r < P r>max . (36) 

Note that this is a non-convex problem. Again, as in the previous subsection, we convert the problem to a 
tractable problem by semi-definite relaxation as follows: 
Problem 5: 

max t 

■R 

s.t. tr((® s (k)-T® N (k))K) >a 2 d r, k = 0, 1, N-l 
tr(* P %) < P r , max 

7e >: 0, (37) 
where $5 (&), &iv(k), and $p(fe) are already defined in Problem 3. 

In Problem 5, the rank constraint rank (1Z) = 1 is dropped by semi-definite relaxation as in Problem 3. Note 
that the relaxed optimization problem is quasi-convex, i.e., for given r the problem is convex. The solution 
of the quasi-convex optimization problem can be obtained by solving its corresponding feasibility problem 
[19]: 

Find 1Z 

s.t. tr -T& N (k))K) > ajr, k = 0, 1, . . . , N - 1 
tr ($ P 7e) < P r , max 

U y 0. (38) 

The feasible set in Problem 5 is convex for any value of r. Let r* be the optimal value of Problem 5. Then, 
we can find the solution to Problem 5 by using the fact that the feasibility problem (38) is feasible for r < r*, 
whereas it is not feasible for r > r*. Based on this, we propose a simple bisection algorithm to solve Problem 
5 as follows: 
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Algorithm 1: 

Step 1: Choose some appropriate interval s.t. r* £ (tl,tr). 
Step 2: Set r = (r L + tr)/2. 

Step 3: Solve the feasibility problem (38) for r. If it is feasible, tl = t. Otherwise, tr = r. 
Step 4: Repeat Steps 1 to 3 until (tr — tl) < e. 

Here, e is the allowed error tolerance for r. Note that the above feasibility problem is a standard SDP problem, 
which can be solved easily by the interior point method [17]. Due to the relaxation, the matrix 1Z obtained 
by solving the relaxed optimization problem may not have rank one in general. In such a case, randomization 
techniques can be applied to find a rank-one solution. 

B. Joint optimization of source power allocation and relay filter 

In the previous subsection, we considered the FF relay design under given source power allocation. How- 
ever, the overall system performance can be improved further by jointly optimizing the source power alloca- 
tion and the FF relay filter. In this subsection, we consider joint source power allocation and FF relay filter 
design to maximize the worst subcanier SNR subject to total source and relay transmit power constraints: 

Problem 6: For given SR channel f, RD channel statistic information (a g ,L g ), FF relay filter order L r , 
maximum available source transmit power P s , ma x, and maximum available relay transmit power P rjmax , 

N-l 

max min SNR k s.t. V] P ak < P s , max and P r < P r . max . (39) 

Ps,0,- ,P B ,N-1,T fcG{0,-,AT-l} f-' 

Note that Problem 6 is a complicated non-convex optimization problem, which does not yield an easy so- 
lution. To circumvent the difficulty in joint optimization, we approach the problem by using alternating 
optimization. That is, first the source power allocation is initialized properly and Problem 4 is solved to op- 
timize the relay filter for given source power allocation. (This step can be done by Algorithm 1.) Then, with 
the given relay filter tap coefficients obtained from Algorithm 1 , the source power allocation is optimized, 
and the alternating procedure is iterated until it converges. Here, the source power allocation problem for a 
given relay filter is formulated as follows: 

Problem 7: For given FF relay filter r, SR channel f , RD channel statistic information (a g , L g ), maximum 
allowed source transmit power P Stmax , and maximum allowed relay transmit power P r ,max, 

N-l 

max min SNR fc s.t. V P Sjk < P s , ma x, P r < ^r.max, and SNR fc > t V k, (40) 

p s ,o,-,p s ,n-i fee{o,-,JV-i} ^ 

where tq is the current maximum worst subcarrier SNR value obtained by Algorithm 1. 
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By introducing a slack variable r and using (30) and (31), Problem 7 can be rewritten as follows: 

max t (41) 

-Ps, Or" ,Ps,N-l,T 



N-l 



s.t. 



^ ^ Ps,k — ^s,max5 (42) 
k=0 

r^EilTEf r < P r , max , (43) 
iVP^r^EiKfcEfr 



> T fe = o, 1, N-l, (44) 

t > r . (45) 

Note that without the relay power constraint (43), the above problem is a simple linear programming (LP) 
with respect to P S; o, • • • , P s ,n-i and r. Indeed, the problem is an LP since the relay power constraint can 
also be written as a linear form in terms of P s ^, as shown below. The relay power (31) can be rewritten as 

E{tr(y t yf)} = tr(RFE{x s xf }¥ H TL H ) + tr(o?RR H ), 
= tr(E{x s xf}F H R^RF) +tr(^RR^), 
( = } tr(Wjv E{ss H } W^F ff R H RF) + tr(cj?RR H ), 

N-l 

( = } Yl P ^ tr ( W^F^R^RF W N ) + tx(a?KR H ) (46) 

k=0 

where is defined in (29), and Wjy is the cyclic prefix extended IDFT matrix given by 

Wat = [wjv-l, Wat„ 2 , • • • , Wo, WjV-1, • • • , WN-Lg-Lr-Lf+3] 1 ■ 

Here, (a) can be verified by using (11) and (b) is due to the assumption of s[k] ~ CJ\f(0, P s ,k) for k = 
0, 1, • • • , N — 1. Using the new expression (46) for the relay transmit power, we obtain an LP optimization 
problem for the source power allocation from the problem (41-45) as 

max r (47) 

Ps,0,— ,Ps,N-l,T 



N-l 

S-t- ^ ^ Ps,k — Ps,m; 



k=0 
N-l 



Y J P s,kCl{k)+C 2 <P r , 



r,max 
k=0 

Ps,kC 3 (k)>T , k = 0, 1, N-l, 
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where C x {k) = tr( e k+1 e T k+1 W^R^RF W N ), C 2 = tr^RR"), and C 3 (k) = ^ggg^ . 
Now, the proposed alternating optimization for joint source power allocation and relay filter design problem 
(39) to maximize the worst subcarrier SNR is summarized as follows. 

Algorithm 2: Given parameters: f, (a g , L g ), L r , P s , max , and P r>max . 
Step 1: Initialize P s>fc for k = 0, • • • , N - 1. 

Step 2: Solve Problem 5 by using Algorithm 1. Set tq as r* obtained from Algorithm 1. 

Step 3: For given r and tq from Step 2, solve Problem (47) to obtain new P s> k, k = 0, • • • , JV — 1. 

Step 4: Go to Step 2 with tl being the solution to Problem (47) in Step 3. 

Step 5: Repeat Steps 2 to 4 until \tq — tl\ < e- 

Here, e is the allowed error tolerance for r. The convergence of the proposed iterative algorithm can easily be 
verified by the monotone convergence theorem applied to t*. Although convergence to the global optimum 
is not guaranteed, it will be seen in the next section that the proposed joint design approach improves the 
system performance significantly. 

IV. Numerical results 
In this section, we provide some numerical results to evaluate the performance of the FF relay design 
methods proposed in Section III. We considered a relay network with an OFDM transmitter, an FF relay, and 
a destination node, as described in Section TL The number of OFDM subcarriers was set as N = 32 with 
minimal cyclic prefix covering the overall FIR channel length in each simulation case. The SR channel length 
and RD channel length were set as Lf = L g = 3, and both SR and RD channel tap coefficients /;'s and #/'s 
were generated i.i.d. according to a Rayleigh distribution, i.e., /; * ~ d ' CN(0, 1) for I = 0, 1, • • • , Lf — 1 and 
gi l ' l rd" Cj\f(0, 1) for I = 0, 1, • • • , L g — 1. Throughout the simulation, the relay and the destination had the 
same noise power of = o\ = 1 and the source transmit power was 20 dB higher than the noise power. Fig. 
3 shows the performance of the first FF relay design method, provided in Problem 3, minimizing the relay 
transmit power subject to the required SNR constraints on subcarrier channels. We chose X = {0, 1, • • • , 27} 
from 32 subcarriers. It is known that for a set of randomly realized propagation channels, it is not easy to 
always guarantee the desired SNR for each subcarrier channel when the desired SNR value is high [7]. Thus, 
in Fig. 3, the line was plotted when Problem 3 was feasible for more than 50 % out of 1000 random channel 
realizations for the given minimum required SNR value for all subcarrier channels in X, and the plotted value 
is the relay transmit power averaged over feasible channel realizations. It is seen that the required relay 
transmit power for the same minimum SNR required by the FF relay is significantly reduced when compared 
with that required by the AF relay. Fig. 4 shows the relay transmit power versus the relay filter length L r for 
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various desired minimum SNR values. It is seen that the required relay transmit power for the same desired 
minimum SNR decreases monotonically with respect to L r , and the FF relay achieves most of the gain only 
with a few taps. Next, we evaluated the performance of the second FF relay design method, provided in 
Algorithm 1, maximizing the worst subcarrier SNR subject to a relay transmit power constraint. Figures 5 
and 6 show the worst subcarrier SNR (averaged over 500 channel realizations) versus the FF relay transmitted 
power for various relay filter lengths and the worst subcarrier SNR (averaged over 500 channel realizations) 
versus the relay filter length L r for various relay transmit power values, respectively. As in the previous case 
of relay transmit power minimization, the gain by the FF relay over the AF relay is significant and most of 
the gain is achieved by only a few filter taps for the FF relay. Recall that our FF relay design is based on 
the channel statistic information for the RD channel. Even exploiting the channel statistic (not the channel 
state) for the RD channel yields significant gain over the AF relay. Now, let us examine the robustness of the 
proposed FF relay design, Algorithm 1, against channel static mismatch. Still, the i.i.d. RD channel model 
with L g = 3 and a l g = \ for all taps was used for Algorithm 1 to obtain the relay filter, but the true channel 
was generated randomly according to a different channel statistic, i.e., a different L g and/or different channel 
power profile. Fig. 7 shows the worst subcarrier SNR (averaged 500 over channel realizations) versus the 
relay transmit power in the case of channel static information mismatch. Here, the worst subcarrier SNR in 
Fig. 7 was computed based on the true channel state and the relay filter obtained by Algorithm 1 with wrong 
channel static information. It is seen in the figure that the proposed FF relay design method is robust against 
the channel static mismatch, and the proposed FF relay outperforms the AF relay even in this case. Finally, 
we investigated the performance of the joint source power allocation and FF relay filter design to maximize 
the worst subcarrier SNR, provided in Algorithm 2, and the result is shown in Fig. 8. It is seen that the joint 
optimization method significantly outperforms the optimization of only the relay filter presented in Algorithm 
1. 

V. Conclusion and Discussion 
In this paper, we have considered the FF relay design for OFDM systems for transparent relay operation 
under two design criteria of minimizing the relay transmit power subject to per-subcarrier SNR constraints 
and maximizing the worst subcarrier SNR subject to a total relay transmit power constraint. We have shown 
that the two design problems reduce to tractable SDP problems and that the FF relay designed by solving 
the obtained SDP problems significantly outperforms the AF relay under both criteria. In closing, we dis- 
cuss several issues and the direction of our future research. First, we assumed channel state information for 
the SR channel and channel statistic information for the RD channel. Although this assumption originates 
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from the consideration of cheap transparent relay operation for FF relays, the assumption of channel statistic 
information for the RD channel makes the proposed FF relay design method useful even for broadcasting pur- 
poses when the destination nodes in the relay cell have similar channel statistics. Second, we assumed SISO 
OFDM systems. However, most current OFDM systems employ MIMO communications, and extension to 
the MIMO case is left as a future work. Third, FF relay design for rate maximization is also an interesting 
problem. This problem can be handled by a similar approach to that used for the rate maximization in single 
carrier systems in [10]. See the end of the Appendix for this issue. 

Appendix 

Proof of Theorem 1 

By introducing a slack variable t, we convert Problem 3 to the following equivalent problem: 

min t (48) 

T,K 

s.t. tr(*pft) < r, (49) 

\x([* s (k)-7 k * N (k)]K) >a 2 dlk , k = 0,1,..., N-l, (50) 

r > 0, (51) 

n y o. (52) 



The Lagrange dual function for the above problem is given by 

N-l 



/ \ 

JV-l N-l 



g{\,{fjL k },v,V) = inf 

T,7c 



k=0 k=0 



V =:Q(A,W) / 

(53) 

where A > 0, > 0}, v > and * ^ are the dual variables associated with (49), (50), (51), and 
(52), respectively. If 1 — A — v ^ or Q(A, — * 7^ 0, then the dual function value is minus infinity 

or we have trivial solutions r = and/or 1Z = 0. Thus, for the nontrivial feasibility of r and 1Z, we have 
1 — A = v (> 0) and Q(A, {fik}) = & (t: 0). Then, the Lagrange dual function is easily obtained as 
g(X, {fJ<k}, v, = J2k=o ^kjk&d an< i the corresponding dual problem is given by 

N-l 



max V fikjko-j 



fc=o 

s.t. 0<A<1, Q(A,{/i fc })>:0, v k >0, k = 0, 1, ...,/Y-l. (54) 

Let A*, {m^}, t*, Q* and 1Z* be the optimal values for the problem, (u* and 1 J r * are automatically deter- 
mined based on these quantities. The dependence of Q on A and is not shown explicitly for notational 
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simplicity from here on.) From the complementary slackness conditions for (49) and (50), we have 

N-l 



A* (tr - r*) + £ {Wd ~ tr ([* s (n) - lk ® N (k)]K*)) = 0, (55) 

which is equivalent to 



A.-0 



<N-1 



£ Vilk<? 2 d - A*r* + tr 



>k=0 



N-l 



{A** P - £ a£[*s(*) - 7fc*iv(fc)]}^ 



fc=0 



0. (56) 



Since the problem (48-52) is a convex optimization problem, the duality gap is zero, i.e., 2^2k=o MfeTfc "! = 
t*. Thus, both the first and second terms in the left-hand side (LHS) of (56) are nonnegative since 2~2k=o t JL k'Yk a d~ 
X*t* = r*(l — A*) > and tr(Q*7£*) > 0. (The trace of the product of two positive semi-definite 
matrices is nonnegative [20].) Therefore, A* = 1 and tr(Q*7£*) = 0. It is obvious that Q* )f for a 
nontrivial TV from tr(Q*7£*) = 0, i.e., the L r x L r matrix Q* does not have full rank. This is because 
tr(Q*7£*) = J2i a it T (Q* VL i u i) = 2~2i a i( u i Q* u i)> where 1Z = J2i °"? u i u ^ lS the eigen-decomposition of 
TV. (If Q* >- 0, then tr(Q*TZ*) > 0.) Note from (56) that 

N-l 

Q* = & P - £ /x£ - 7fc*tf(*0) , (57) 

fc=o 

where <&p is a positive definite matrix defined in (31), «J?at(A;) is a positive semi-definite matrix defined in 
(28), and defined in (26) is a rank-one matrix by Lemma 2. Now, under the assumption that all the 

SNR constraints except one are satisfied with strict inequality, we have /Zj 7^ for some i and fij = 0, Vj 7^ i 
from the complementary slackness conditions. In this case, Q* is given by 

Q* = <& P + rift&NQ - tfGsii) ■ (58) 
rank L r rank 1 

Due to the structure of Q* in (58), the rank of Q* is larger than or equal to L r — 1. Since Q* )/- 0, 
rank(Q*) = L r - 1. Since = tr(Q*ft*) = Efef 1 »»( v f ^* v <) ( where Q* = Ef=i"SiVivf is the 
eigen-decomposition of Q*), we conclude that 1Z* has nullity L r — 1 and thus has rank one under the as- 
sumption of Theorem 1 . ■ 

Lemma 2: If N — Lt + 1 > Lt + L g — 1, <&s(fc) has rank one regardless of the value of k. 
Proof of Lemma 2: Recall that (see (26) and (33)) 

(NP^y^sik) = EiKfeEf , K k = I Lg ®K k , K k = F*T*w fc wf T T F T (59) 
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Let Ei in (27) be partitioned as Ei = E^ 1} , . . . , e[ N+L "~^ . Note that F and T are Toeplitz 

matrices and that w^wF is also a Toeplitz matrix regardless of k due to the property of DFT matrices. It is 
not difficult to show that K^(l : N — Lf + 1, 1 : N — L* + 1) is a Toeplitz matrix, where A(a : b,c : d) 
denotes a submatrix of A composed of the rows from a to b and columns from c to d. Now, can be 

rewritten as 

(JVP^oJ)- 1 *^*) = E«K fc E« H + E< 2) K fc E? )fl + ■■■ + n[ Ls) K k B[ L9)H ■ (60) 

Here, the operation E^K^E}^ extracts a L r x L r submatrix K/-(i : L r + i — 1, % : L r + i — 1) from 
Kfe. If iV — L r + 1 > L r + L g — 1, this operation extracts the same submatrix from regardless of i since 
Kfc(l : N — L r + 1, 1 : N — L r + 1) is a Toeplitz matrix. Thus, we have 

(JVP^r^sCfc) = L 9 ES 1} K fc ES 1)H = L 9 (ES 1) F*T*w fc )(Ef ) F*T*w fc ) // , (61) 

and &s(k) nas ran k one if the condition N — Lf + 1 > L t + L g — 1. ■ 

7omf Optimization of Source Power Allocation and FF Relay Filter for Rate Maximization 

This problem can be formulated by using the results in this paper, and the problem is expressed as follows: 
Problem 8: For given f, (L g) a 2 ), L r , P S:max and P r ,max, 



max 

fa, Or" ,Ps,N- 



.... S2 b ^ 1+ ^ r «E 2 M fc E fr + ^j <62) 



AT-1 



s.t. < P s , max , (63) 

fc=0 

r H Ein*Ef r < P r , max . (64) 

The problem (62-64) is not a convex problem and thus, does not yield a closed-form solution. However, this 
problem is a discrete frequency version of the problem considered in [10]. Thus, the approach in [10] based on 
the projected gradient method can be applied with the new design variable u := [P s> o,--- ,P s ,N-i,fo,ri,--- ,ri r 
and (p(u) := Ylk=o 1 1°§ ( 1 + N „21h ^ ^ ~^?~ri~2~ ) • This a PP roacn yields a nontrivial rate gain over the 



AF scheme. Please see [10] and [21] for detail. 
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Fig. 1 

The considered relay network 
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Fig. 2 

The considered relay network with an OFDM transmitter, an FF relay, and a destination node 
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